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We investigate several Pb(n, n'7) and Ge(n, n'7) reactions. We measure 7-ray production from 
Pb(n, n'7) reactions that can be a significant background for double-beta decay experiments which 
use lead as a massive inner shield. Particularly worrisome for Ge-based double-beta decay experi- 
ments are the 2041-keV and 3062-keV 7 rays produced via Pb(n, n'7). The former is very close to 
the 76 Ge double-beta decay endpoint energy and the latter has a double escape peak energy near 
the endpoint. We discuss the implications of these 7 rays on past and future double-beta decay 
experiments and estimate the cross section to excite the level that produces the 3062-keV 7 ray. 
Excitation 7-ray lines from Ge(n, n'7) reactions are also observed. We consider the contribution of 
such backgrounds and their impact on the sensitivity of next-generation searches for neutrinoless 
double-beta decay using enriched germanium detectors. 

PACS numbers: 23.40.-s, 25.40.Fq, 29.40.Wk 
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I. INTRODUCTION 



Neutrinoless double-beta decay plays a key role in 
understanding the neutrino's absolute mass scale and 
particle-antiparticle nature [H, H, [3, 4]. If this nuclear de- 
cay process exists, one would observe a mono-energetic 
line originating from a material containing an isotope 
subject to this decay mode. One such isotope that may 
undergo this decay is 76 Ge. Germanium-diode detectors 
fabricated from material enriched in 76 Ge have estab- 
lished the best half-life limits and the most restrictive 
constraints on the effective Majorana mass for the neu- 
trino [1, § . One analysis Q of the data in Ref. @ claims 
evidence for the decay with a half-life of 1.2 x 10 25 y. 
Planned Ge-based double beta decay experiments [1, Q 
will test this claim. Eventually, these future experiments 
target a sensitivity of > 10 27 y or ~ 1 event /ton- year 
to explore mass values near that indicated by the atmo- 
spheric neutrino oscillation results. 

The key to these experiments lies in the ability to re- 
duce intrinsic radioactive background to unprecedented 
levels and to adequately shield the detectors from ex- 
ternal sources of radioactivity. Previous experiments' 
limiting backgrounds have been trace levels of natural 
decay chain isotopes within the detector and shielding 
components. The 7-ray emissions from these isotopes 
can deposit energy in the Ge detectors producing a con- 
tinuum, which may overwhelm the potential ncutrino- 
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less double-beta-decay signal peak at 2039 keV. Great 
progress has been made identifying the location and ori- 
gin of this contamination, and future efforts will substan- 
tially reduce this contribution to the background. The 
background level goal of 1 cvent/ton-year, however, is an 
ambitious factor of 400 improvement over the currently 
best achieved background level Q. If the efforts to re- 
duce the natural decay chain isotopes are successful, pre- 
viously unimportant components of the background must 
be understood and eliminated. The potential for neutron 
reactions to be one of these background components is 
the focus of this paper. The work of Mei and Himefloj 
recognized that (n, n'j) reactions will become important 
for ton-scale double-beta decay experiments. Specifically, 
we have studied neutron reactions in Pb and Ge, materi- 
als that play important roles in the Majorana |8[ design. 
But since lead is used by numerous low-background ex- 
periments, the results will have wider utility. 

This paper presents measurements and simulations of 
Pb(n, n'"f) and Ge(n, n'7) reactions and estimates the re- 
sulting background for Ge-detector based, double-beta 
decay experiments for a given neutron flux. With these 
results, we then use the neutron flux, energy spectrum, 
angular distribution, multiplicity and lateral distribu- 
tions determined in (Io| to estimate the background in 
Ge detectors situated in underground laboratories. In 
Section[TT]we describe the experiments, the data, and the 
simulations. In Scctions llIII and llVl we describe the analy- 
sis of these data. Section [TVl also discusses the important 
Pb(n, n'7) production of 7 rays at 2041 and 3062 keV. 
The former is dangerously near the 2039-keV Q-value for 
zero- neutrino double-beta decay in 76 Ge and the latter 
can produce a double-escape peak line at 2040 keV. These 
dangerous processes for Ge-based double-beta decay ex- 
periments are discussed for the first time in this work. 
Section [V] determines an overall background model for 
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our detector and the implications of this model for fu- 
ture experimental designs. It also considers the relevant 
merits of Cu versus Pb as shielding materials, and the 
use of depth to mitigate these backgrounds is discussed. 
We also consider the possibility that the double-escape 
peak of the 3062-keV 7 ray could contribute to the signal 
claimed in Ref. Q- Finally, we summarize our conclu- 
sions in Section EH 



II. THE MEASUREMENTS 

We collected five data sets to explore the implications 
of (n, n'j) for double-beta decay experiments. All mea- 
surements were done in our basement laboratory at Los 
Alamos National Laboratory. The laboratory building is 
at an atmospheric depth of 792 g/cm 2 and provides about 
1 mwc concrete (77 g/cm 2 ) overburden against cosmic 
ray muons. 

Three data sets were taken with a CLOVER detec- 
tor This detector is a set of 4 n-type, segmented 
germanium detectors. The four crystals have a total nat- 
ural germanium mass of 3 kg and each crystal is seg- 
mented in half. The CLOVER detector and its opera- 
tion in our laboratory were described in Ref. [l^ . The 
remainin g tw o measurements were done with a PopTop 
detector |13( set up in coincidence with a Nal detector. 
The PopTop is a 71.8-mm long by 64- mm diameter p- 
type Ge detector. Taking into account the central bore, 
the detector is 215 cm 3 or 1.14 kg. The Nal crystal is 
15.25-cm long by 15.25-cm diameter and is directly con- 
nected to a photo-tubc. All data were read out using 
a pair of X Ray Instrumentation Associates (XIA) [14[ 
Digital Gamma Finder Four Channel (DGF4C) CAMAC 
modules. The CAMAC crate is connected to the PCI bus 
of a Dell Optiplcx computer running Windows 2000. The 
system was controlled using the standard software sup- 
plied by XIA. This data acquisition software runs in the 
IGOR Pro environment [la ] and produces binary data 
files that were read in and analyzed using the ROOT 
framework [l6j. 

The data sets include: 

1. A background run with the CLOVER 

2. A Th-wirc source run with the CLOVER 

3. An AmBe source run with the CLOVER using two 
different geometries of moderator 

4. An AmBe source run with the PopTop surrounded 
by lead 

5. An AmBe source run with the PopTop surrounded 
by copper 

In this section we describe the experiments and the data 
collected. 



A. The Experimental Configurations 

The CLOVER was surrounded by 10 cm of lead shield- 
ing to reduce the signal from ambient radioactivity. Un- 
derneath and above the lead was 5 cm of 30%-loaded 
borated polyethylene to reduce thermal neutrons. The 
background run done in this configuration lasted 27.13 
live-days. The configuration for the Th-source run was 
similar, but with some lead removed to expose the detec- 
tor to the source. The Th source run had a live time of 
1337 seconds. 

The setup was modified somewhat from this 
background-run configuration for the measurements with 
the AmBe source. Fig. []] shows the configuration for 
one of the AmBe measurements. For these data, the 
CLOVER was shielded on four sides with 10 cm of lead. 
The AmBe source, 30 mCi of 241 Am with a calibrated 
neutron yield of ss 63,000 Hz (±0.7%), was on one side 
of the CLOVER with 5 cm of lead and a layer of pure 
polyethylene moderator (either 10 or 15 cm thick) be- 
tween the source and detector. The data acquisition sys- 
tem is inactive during data transfer. Only the AmBe 
runs had a large enough event rate for the dead time to 
be appreciable. A 6.13-h live-time data run (57% live) 
was taken with 15 cm of moderator, and another 3.57-h 
live-time data (38% live) run was taken with 10 cm mod- 
erator (pictured). For the analysis presented below, the 
data from these two configurations were combined, and 
thus the AmBe-CLOVER data set contains 9.7 h of live 
time. The observed energy spectrum extended from sa 
10-3100 keV for these data sets. 

During the analysis of the AmBe data, we observed 
a weak line at 3062 keV. This energy corresponds to a 
7-ray transition in 207 Pb, and we therefore hypothesized 
that it was generated via Pb(n, n'7). The double-cscapc- 
peak (DEP) energy (2040 keV) associated with this 7 ray 
is very dangerous for 76 Ge neutrinoless double-beta de- 
cay experiments because it falls so close to the transition 
energy (2039 keV). Furthermore because the DEP is a 
single-site energy deposition, it cannot be distinguished 
from double-beta decay through event topology. This is 
in contrast to a full-energy 7-ray peak, which tends to 
consist of several interactions and therefore is a multiple- 
site deposition. (See [l2| for a discussion of the use of 
event topology to reduce background in Ge detectors.) 

The final two measurements were intended to study 
this 3062-keV line in the spectrum and demonstrate its 
origin. In both cases a PopTop Ge detector faced a 15.25 
cm by 15.25 cm Nal detector for coincidence data. By se- 
quentially placing a Pb and then a Cu absorber between 
an AmBe source and a PopTop Gc detector, we tested 
the hypothesis that the line was due to neutron interac- 
tions in Pb. By looking for an coincident energy deposit 
in the Nal detector, we could be assured the Ge detector 
signal originated from a neutron interaction in the sam- 
ple. An energy deposit threshold in the Nal of greater 
than 200 keV was required for a coincidence. The Pop- 
Top was placed 27.3 cm from the Nal detector with the 
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source placed 20.3 cm (7 cm) from the Ge (Nal) detector. 
For the lead study, 5 cm of lead was placed directly be- 
tween the Ge detector and the source. Additional lead, in 
the form of 5-cm-thick bricks was positioned around the 
4 sides of the Ge detector to reduce room background. 
For the copper study, a 0.5-cm thick Cu tube was placed 
around the PopTop and a 5-cm Cu block was placed be- 
tween the PopTop and the source. For this final run, all 
the lead was removed. For both of these sets of data, the 
observed spectra extended from « 125 keV to « 9 McV. 
For the PopTop data, the Pb and Cu runs were of 19.12 
h and 17.76 h live-time, respectively. 



spectra are due to the sources and other contributions 
can be safely ignored. For example, the 2614.5-keV line 




FIG. 1: The CLOVER detector as configured for the AmBe 
source run. The setup at the time of this photograph used 
4" of polyethylene. One wall of the lead shield was removed 
only to clarify the relationship between the AmBe source, 
moderator, and the CLOVER. 



B. The Data sets 

The crystals were individually calibrated, and the re- 
sulting spectra summed together to form a single his- 
togram. The peaks within each of the 3 CLOVER data 
sets were identified and their intensities determined. If an 
event had 2 crystals that responded in coincidence, the 
histogram would have two entries. Therefore the spectra 
we analyzed and simulated included all single-crystal en- 
ergy deposits. By not eliminating events that registered 
signals in more than 1 of the CLOVER Ge detectors, 
we maximized the event rate. The peak strengths were 
estimated by fitting a Gaussian shape to peaks and a 
flat background to the spectrum in the region near the 
peak. For the nuclear recoil lines, the peak shape was 
assumed to be a triangle and not Gaussian. In Table U 
the uncertainties derive from this fit. A summary of the 
peak strengths is given in Table [J and the spectra them- 
selves are shown in Fig. [2j The data sets were chosen 
to help decouple line blcndings. Because the rates in 
all peaks and continua are much higher for the source- 
induced data than for the background, features in those 



can arise from either the decay of 2U8 T1 or 0S Pb(n, n'j). 
When exposed to a Th source, Tl decay dominates the 
spectrum, whereas when exposed to an AmBe source, 
(n, n'7) dominates. Hence by normalizing the rate in 
this line to the rate in a pure neutron-induced transition 
(e.g. the 596-keV 74 Ge(n,n'7)), we can determine the 
relative contribution of the two processes to the back- 
ground spectrum. In fact, in the background data, both 
processes contribute to this line. 

Some comments on our choices for line identification 
are in order. For an isotope such as 72 Ge where a neutron 
capture leads to a stable nucleus, almost all (n,7) lines 
could also be interpreted as (n, 77/7) lines in the result- 
ing nucleus; in this case 73 Ge. For isotopes within the 
detector however, such as the 53.5-keV 72 Ge(n,7) tran- 
sition, the competing 73 Ge(77, 77/7) line would be a sum 
of this 7-ray energy and the recoil nucleus energy. At 
these low energies where the recoil is a fair fraction of 
the 7-ray energy, the (n, 77/7) would simply contribute to 
the continuum and not be observed as a line. For the 
high energy cases, the blend of a mono-energetic 7-ray 
line and a (77,77/7) process might be present. 

For the calibration runs, our threshold was approxi- 
mately 70 keV. Also note that we used a thorium wire as 
a calibration source. Since the wire is pure natural Th, 
we observe the Th X rays in that data. In contrast, the 
background run shows lines from the thorium chain as a 
contaminant, therefore those lines are absent. 

In all spectra, there are a few lines we have not iden- 
tified. 

TABLE I: A summary of the observed lines in the various spectra taken 
with the CLOVER detector. Blank entries indicate that no significant 
peak feature above the continuum was found. Single and double escape 
peaks are labeled by SEP and DEP respectively. Line assignments for 
which we are unsure are indicated by a question mark. Line energies are 
taken from the Table of Isotopes [13] ■ 



Energy 


Process 


Count Rates 


(keV) 




backgrnd 
(per hr) 


Thorium 
(Hz) 


CLOVER 
AmBe (Hz) 


23.4 


' u Ge(n, 7 ) 






1.017(5) 


46.5 


2io pb 


112.75(42) 






53.2 ] 
53.5 


> 


234 U \ 
72 Ge(n, 7 ) / 


61.01(31) 




2.079(8) 


63.2 


234 Th 


93.02(38) 






67.7 


230^ 


24.05(19) 




0.591(4) 


68.8 


72 Ge(n, 7 ) 






0.440(4) 


72.80 


Pb x-ray 




10.5(1) 


0.506(4) 


74.97 


Pb x-ray 


663.5(1.0) 


34.1(2) 


1.703(7) 


76.7 


Unidentified 




29.0(2) 




84.4 ' 




228 Th >, 






84.5 




Pb x-ray > 


115.28(42) 


10.9(1) 


0.930(5) 


84.9 




Pb x-ray J 








87.2 A 












87.4 




Pb x-ray 


55.65(29) 


16.8(1) 


0.262(3) 


87.9 













Continued 
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TABLE I - continued TABLE I - continued 



Energy 


Process 


Count Rates 


Energy 


Process 


Count Rates 


(keV) 




backgrnd 
(per hr) 


Thorium 
(Hz) 


CLOVER 
AmBe (Hz) 


(keV) 




backgrnd 
(per hr) 


Thorium 
(Hz) 


CLOVER 
AmBe (Hz 


89.9 


Th x-ray 




32.4(2) 




478.6 


228 Ac 




0.39(2) 




92.7 


234 Th 


171.38(51) 




0.060(1) 






' 10 B(n,a) 7 Li* 










93.4 
96.0 


Th x-ray 
115 In(n, 7 ) ? 




47.5(2) 


0.166(2) 


470-485 


< 


7 Li*( 7 ) 7 Li 
Doppler 


! 






signf. 


99.5 


228 Ac 


13.71(15) 


3.0(1) 








Broadened 










105.3 


Unidentified 


8.99(12) 






492.9 


73 Ge(n, 7 ) 






0.123(2) 


104.8 
105.6 J 


\ 


Th x-ray 




20.6(1) 




499.9 
503.9 


70 Ge(n, 7 ) 
228 Ac 




0.34(2) 


0.453(4) 


108.7 


Th x-ray 




7.6(1) 




509.3 


228 Ac ) 










109.9 
129.1 


19 F(n,n' 7 ) 

228 * 

Ac 


43.00(26) 
12.91(14) 


3.2(1) 


0.506(4) 


510.7 
510.7 


208rpj 

208 Pb(n,n' 7 ) 


> 


171.93(51) 


16.8(1) 


3.409(10) 


139.7 


74 Ge(n, 7 ) 


47.20(27) 




2.339(8) 


510.9 


Annih. 7 










143.9 


230 Th 


20.03(18) 






516.2 


35 Cl(n, 7 ) 






0.160(2) 


154.0 


228 Ac 


7.69(11) 


1.5(1) 




537.5 


206 Pb(n,n' 7 ) 


5.12(9) 




0.158(2) 


159.7 


77m Qe 






0.114(2) 


562.9 1 


228 Ac \ 




12.83(14) 


1.52(3) 


0.244(3) 


162.4 


115 In(n, 7 ) 


10.71(13) 




1.073(6) 


563.0 J 


76 Ge(n,n' 7 ) j 




174.9 


70 Ge(n, 7 ) 


7.45(11) 




0.763(5) 


569.7 


207 Pb(n,n' 7 ) 


14.17(15) 




0.422(4) 


186.1 
186.2 


■ 


22b Ra 1 
115 In(n, 7 ) j 


> 


114.60(42) 




0.323(3) 


572.3 
574.7 


228 Ac 
74 Ge(n, 7 ) 




0.53(2) 


0.091(2) 


197.1 < 
198.4 




19 F(n,n' 7 ) 
71 Ge sum 


} 


81.04(35) 




2.328(8) 


583.1 


\ 


208 rpj 

208 Pb(n, n'y) 


} 


71.64(33) 


49.6(2) 


0.256(3) 


199.2 
202.6 


228 Ac 
115 In(n, 7 ) 




0.66(2) 


0.061(1) 


595.9 




) 74 Ge(n,n' 7 ) ] 
1 73 Ge(n, 7 ) 




59.90(30) 




1.869(7) 


209.5 


228 Ac 


19.38(17) 


10.9(1) 




608.3 


73 Ge(n, 7 ) 






0.333(3) 


215.5 


228^ 


2.43(6) 


0.92(3) 




609.2 


2i4 Bi 


60.11(30) 






238.6 


212p b 


295.77(67) 


139.9(1) 


0.105(2) 


629.6 


72 Ge{n,n'j) 




0.078(2) 


242.0 


214p b 


57.49(30) 


9.2(1) 




648.2 


115 In(n, 7 ) 






0.025(1) 


247.1 


70 Ge(n, 7 ) 






0.070(1) 


657.2 


206 Pb(n,n' 7 ) 






0.047( 1) 


253.7 


74 Ge(n, 7 ) 


2.76(7) 




0.410(3) 


662.0 


137 Cs 


9.04(12) 






270.2 


228 Ac 


21.39(18) 


9.1(1) 




663.8 


206 Pb(n, n'j) 




0.069(1) 


273.0 


115 In(n, 7 ) 






0.055(1) 


669.0 


70 Ge(n,n' 1 ) 






0.030(1) 


277.4 


f 208 Tl \ 
\ 208 Pb(n,n' 7 ) J 


12.98(14) 


5.4(1) 


0.086(2) 


692.4 
701.0 


72 Ge(n,n'e~) 
74 Ge(n,n' 7 ) 


87.70(37) 




2.406(8) 
0.082(2) 


284.6 


Unidentified 


2.79(7) 






708.2 


70 Ge(n, 7 ) 






0.176(2) 


288.1 


212 Bi ? 




1.03(3) 




727.3 


212 Bi 


15.72(16) 


11.7(1) 




295.2 


214p b 


58.02(30) 






747.7 


70 Ge(n, 7 ) 


0.047(1) 


297.2 


\ 


"Ge(n, 7 ) 1 


► 






0.068(1) 


755.3 


228 Ac 


2.19(6) 


1.52(3) 




298.7 j 




115 In(n, 7 ) j 








763.1 


208 pj 




0.85(3) 




300.1 


2i2 Pb 


18.59(17) 


9.7(1) 




763.1 


208 Pb(n,n' 7 )? 




0.032(1) 


306.2 


70 Ge(n, 7 ) 


1.12(4) 




0.046(1) 


766.6 


224m p a V 


4.85(9) 




321.4 


228 Ac 




0.67(2) 




768.4 


214 Bi ( 






326.0 


70 ' 72 Ge(n,n' 7 ) 






0.487(4) 


771.8 


228 Ac 


2.02(6) 


2.11(4) 




328.3 


228 Ac 


10.02(12) 


8.2(1) 




782.0 


228 Ac 


0.67(2) 




332.9 


228 Ac 




1.14(3) 




785.5 


212 Bi 


4.19(8) 


1.57(3) 




335.5 


115 In(n, 7 ) 






0.028(1) 


786.3 


35 Cl(n, 7 ) 1 








0.041(1) 


338.7 


228 Ac 


45.13(26) 


31.5(2) 




786.8 


208 Pb(n,n' 7 ) 


► 






351.9 


214p b 


95.11(38) 






788.4 


35 Cl(n, 7 ) 1 






0.064(1) 


354.1 


Unidentified 






0.043(1) 


788.7 


7 °Ge(n, 7 ) / 






385.1 


115 In(n, 7 ) 






0.048(1) 


795.0 


228 Ac 


7.92(11) 


6.1(1) 




391.3 


70 Ge(n, 7 ) 






0.053(1) 


798.0 


208 Pb(n, n'-y) 


0.023(1) 


409.8 


228 Ac 


4.05(8) 


4.3(1) 




803.1 


206 Pb(n,n' 7 ) 


20.90(18) 




0.850(5) 


416.9 




2.21(6) 




0.359(3) 


806.2 


2i4 Bi 


3.03(7) 




438.9 


Unidentified 


2.54(6) 






808.2 


70 Ge(n, 7 ) 




0.048(1) 


445.2 


74 Ge(n, 7 ) 






0.037(1) 


818.6 


116m In 






0.064(1) 


452.3 


212 Bi? 




0.82(2) 




824.9 




1.03(4) 




463.3 


228 Ac 


10.96(13) 


9.2(1) 




830.4 


228 Ac 


0.70(2) 




474.0 


72 Ge(n, 7 ) ? 


2.54(6) 






834.1 


72 Ge(n,n' 7 ) 


45.15(26) 


0.290(3) 



Continued Continued 
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TABLE I - continued TABLE I - continued 



Energy 


Process 


Count Rates 


Energy 


Process 


Count Rates 


(keV) 




backgrnd 
(per hr) 


Thorium 
(Hz) 


CLOVER 
AmBe (Hz) 


(keV) 




backgrnd 
(per hr) 


Thorium 
(Hz) 


CLOVER 
AmBe (Hz) 


835.6 


228 Ac 




2.30(4) 




1261.0 


Y4 Ge(n, n' 7 ) 






0.019(1) 


840.4 


228 Ac 




1.19(3) 




1281.0 


2i4 Bi 


0.74(3) 






843.8 


27 Al(n,n' 7 ) 


4.48(8) 




0.112(2) 


1286-7 


228 Ac Blend 




0.14(1) 




846.9 


76 Ge(rj,n' 7 ) 






0.062(1) 


1293.5 


116m In 


4.09(8) 




0.462(4) 


860.4 
860.4 


Tl 1 

208 Pb(n,n' 7 ) J 




8.64(12) 


6.0(1) 


0.090(2) 


1298.8 
1332.5 


70 Ge(n, 7 ) 
74 Ge(n,n' 7 ) 






0.087(2) 
0.018(1) 


865.0 


Unidentified 






0.094(2) 


1332.5 


60 Co 


12.27(14) 






867.9 


73 Ge(n, 7 ) 


4.25(8) 




0.466(4) 


1344.5 




74 Ge(n, 7 ) 










881.0 


206 Pb(n, n' 7 ) 


2.50(6) 




0.151(2) 


1345.9 


1 


206 Pb(n,n' 7 ) 


1 


0.52(3) 




0.017(1) 


892.9 


2i2 Bi 




0.42(2) 




1347.7 




7 °Ge(n, 7 ) 










894.3 


72 Ge(n,n' 7 ) 






0.029(1) 






f 228 Ac sum \ 










897.8 


207 Pb(n, n' 7 ) 


6.28(10) 




0.199(2) 


1374.2 




I 964 + 409 \ 






0.28(1) 




904.1 


228 Ac 




0.94(3) 








1 911 +463 1 










911.2 


228 Ac 


48.86(27) 


37.1(2) 




1378.0 


Bi 


3.24(7) 






934.1 


214 Bi 


1.99(6) 






1378.8 


,u Ge(n, 7 ) 






0.065(1) 


958.4 


228 Ac 




0.37(2) 




1393.8 


206 m-./ / \ 

Pb(n, n 7 ) 






0.016(1) 


960.9 


74 Ge(n,n' 7 ) 






0.095(2) 


1401.5 


214 Bi 


0.58(3) 






964.4 


228 Ac 


11.22(13) 


6.2(1) 




1408.6 


214 Bi 


1.84(5) 






968.8 


228 Ac 


26.83(20) 


21.6(1) 




1413.6 


73 Ge(n,n' 7 ) 






0.018(1) 


981.0 


206 ' 8 Pb(n,n' 7 ) 






0.035(1) 


1431.1 


228 Ac 




0.15(1) 




988.4 


228 Ac 


1.95(5) 


0.20(1) 




1433.5 


206 Pb(n, n'i) 






0.020(1) 


993.7 


\ 


' 74 Ge(n,n' 7 ) 


I 






0.026(1) 


1436.9 


208 Pb(n,n' 7 ) 






0.017(1) 


995.1 




206 Pb(n, n'j) 








1459.2 


228 Ac 




0.80(2) 




999.5 


74 Ge(n,n' 7 ) 






0.034(1) 


1461.0 


40 K 


30.18(22) 




0.066(1) 


1001.5 


224m p a 


8.03(11) 






1463.9 


72 Ge(n,n' 7 ) 






0.114(2) 


1004.5 


228 Ac 




0.17(1) 




1466.8 


206 Pb(n,n' 7 ) 






0.032(1) 


1014.5 


27 Al(n,?i' 7 ) 


7.46(11) 




0.173(2) 


1471.6 


73 Ge(n, 7 ) 






0.047(1) 


1033.1 


228 Ac 




0.20(1) 




1489.2 


74 Ge(n,n' 7 ) 






0.025(1) 


1040.1 1 
1040.8 J 


70 Ge{n,n'j) 


16.89(16) 




0.210(2) 


1496.2 
1501.7 


228 Ac 
228 Ac 


0.73(3) 


0.85(3) 
0.42(2) 




1063.7 


207 Pb(n,n' 7 ) 


8.47(11) 




0.145(2) 


1508.9 


116m In 






0.069(1) 


1065.0 


228 Ac 




0.47(2) 




1508.9 


2i4 Bi 


1.95(3) 






1078.8 


212 Bi 


1.11(4) 


0.62(2) 




1512.7 


212 Bi 




0.38(2) 




1093.9 


208 T1 sum 511+583 




0.90(3) 




1538 


214 Bi 


0.54(3) 






1095 




' 207 Pb(n,n' 7 ) 










1557.1 


228 Ac 




0.14(1) 




1095.8 


! 


™Ge(n, 7 ) 


) 


7.81(11) 




0.464(4) 


1580.8 


228 Ac 


0.68(3) 


0.55(2) 




1096.9 




116m Tn 










1588.3 


228 Ac 


6.06(10) 


2.94(5) 




1101.3 


74 Ge(n,n' 7 ) 






0.123(2) 


1592.5 


208 T1DEP ] 










1105.6 


74 Ge(n,n' 7 ) 






0.019(1) 


1592.5 


208 Pb DEP 




7.24(11) 


2.10(4) 


0.107(2) 


1110.4 


228 Ac sum 




0.52(2) 




1593.0 


207 Pb(n,n' 7 ) 










1120.6 


214 Bi 


11.70(13) 






1599.3 


214 Bi 


0.38(2) 






1122.5 


228 Ac sum 




0.26(1) 




1601.1 




' 35 Cl(n, 7 ) ~ 








0.018(1) 


1126 


i 


208 Pb(n,n' 7 ) 


} 






0.022(1) 


1602.0 


\ 


74 Ge(n,n' 7 ) 


\ 








^ 72 Ge(n,n' 7 ) 








1614.9 


208 Pb(n,n' 7 )? 






0.020(1) 


1131.6 


73 Ge(n, 7 )? 


1.34(5) 




0.034(1) 


1620.5 


2i2 Bi 


1.81(5) 


1.32(3) 




1139.4 


70 Ge(n, 7 ) 


1.34(5) 




0.077(2) 


1625.0 


228 Ac 




0.34(2) 




1153.5 


228 Ac 




0.15(1) 




1630.7 


228 Ac 


1.82(5) 


1.46(3) 




1155.2 


2i4 Bi 


1.08(4) 






1631.5 


74 Ge(n,n' 7 ) \ 








0.021(1) 


1164.9 




f J5 Cl(n,n' 7 ) 1 




0.49(3) 




0.072(1) 


1632.0 


70 Ge(n, 7 ) / 








1166.0 




1 72 Ge(n,n' 7 ) j 


> 




1634.0 


76 Ge(n, 7 ) 






0.016(1) 


1173.5 


60 Co 


16.01(14:) 






1638.3 


228 Ac 


0.43(3) 


0.41(2) 




1201.2 


p(n, 7 )d DEP 






0.415(3) 


1640.4 


208 Pb(n,n' 7 ) 








0.026(1) 


1204.2 




( "Ge(n, 7 ) 1 
1 74 Ge(n,n' 7 ) 1 


* 


9.39(12) 




0.163(2) 


1640 
1661.3 


74 ' 76 Ge(n,n' 7 ) 
214 Bi 


} 


0.36(2) 




1226.7 


74 Ge(n,n' 7 ) 






0.017(1) 


1666.3 


228 Ac 




0.17(1) 




1238.4 


214 Bi 


5.22(9) 






1699.5 


206 Pb(n,n' 7 )? 






0.021(1) 


1246.9 


228 Ac 




0.58(2) 




1704.5 


206 Pb(n,7i' 7 ) 


0.72(3) 




0.041(1) 



Continued Continued 
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TABLE I - continued 



Energy 


Process 


Count Rates 


(keV) 




backgrnd 


Thorium 


CLOVER 






(per hr) 


(Hz) 


AmBe (Hz) 


1710.9 


™Ge(n, n'7) 






0.327(3) 


1712.2 


p(n,7)d SEP 






0.156(2) 


1725.7 


207 Pb(n, 717) 






0.029(1) 


1729.6 




3.59(7) 






1764.7 


21 4n- 


-1-1 /^o/io\ 

11.68(13) 






1779.0 


27 Al(n/y) 8 A1 \ 
28 Al =^ 28 Si J 


2.13(6) 




0.127(2) 


1806.0 


212 Bi 




0.11(1) 




1844.5 


2Ub Pb(n,n 7) 






0.044(1) 


1846.9 


2i4 Bi 


2.26(6) 






1940.4 


74 Ge(rc,n' 7 ) 






0.027(1) 


1951.1 


35 Cl(n, 7 ) 






0.032(1) 


1959.3 


35 Cl(n, 7 ) 






0.023(1) 


2092.1 


206 Pb(n, n'j) \ 






0.039(1) 


2092.7 


207 Pb(n,n'7) J 






2103.8 


/ 208 Tl SEP T 
\ 208 Pb SEP J 


5.21(9) 


2.25(4) 


0.090(2) 


2112.1 


116 "Tn 






0.061(1) 


2118.5 


214 Bi 


0.64(3) 






2204.0 


214 Bi 


3.73(8) 






2223.3 


p(n,7)d 






5.813(13) 


2390.5 


116m In 






0.015(1) 


2448.5 


2i4 Bi 


0.51(3) 








f 208 Tl 1 


16.3(1) 


0.729(5) 


2614.5 


\ 208 Pb(n,n' 7 ) / 


39.39(25) 


2650.3 


206 Pb(n,n' 7 )? 






0.011(1) 


2686 


sum 208 Tl? 




0.12(1) 




2892 


sum 208 Tl 




0.08(1) 




3061.9 


207 Pb(n,n'7) 






0.010(1) 



C. Neutron Spectra Simulation 

Fast neutrons (from 100 MeV to 1 GeV or more) tend 
to produce additional neutrons through nuclear reactions 
as they traverse high-Z material. In particular the flux 
of neutrons will increase several-fold while the average 
neutron energy decreases through these processes. As a 
result, fast neutrons will penetrate deep into a shield pro- 
ducing additional neutrons at lower energies. These low 
energy neutrons (~< 20 MeV) give rise to a substantial 
7-ray flux because (n,n 7) cross sections are large near 
10 MeV, but become small at higher energies. Hence 
it is these secondary lower-energy neutrons that inter- 
act with the shield and detector materials to produce 7 
rays, which can give rise to background in double-beta 
decay experiments. To understand the process by which 
high energy neutrons influence the low-energy neutron 
flux and, in turn, the observed 7-ray flux, we simulated 
neutrons impinging on an outer shield and tracked how 
their spectrum changed as the particles traversed the 
shield. We also simulated the production of neutron- 
induced 7 rays and how the Ge detector responded to 
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FIG. 2: The AmBe and background spectra taken with the 
CLOVER. 



them. Specifically, wc performed simulations of several 
geometries including: 

1. A simulation of the cosmic-ray produced neutrons 
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with energy up to 1 GeV at our lab in Los Alamos 
and their propagation through a 10-cm Pb shield. 
The response of the CLOVER detector to 7 rays 
produced by neutron interactions in the shield was 
simulated. This simulation, compared to our back- 
ground data, tests the precision to which we can 
model neutron production, scattering with sec- 
ondary neutron production and (n,n'7) interac- 
tions. 

2. A simulation of the neutron flux induced on the 
CLOVER from the AmBe source (neutron energy 
<11.2 MeV) passing through 15 cm of polyethylene 
before impinging on the 10-cm Pb shield. Since the 
neutron flux is of low energy, this simulation tests 
the precision to which we model (n,n'7) interac- 
tions. 

3. A simulation of the neutron flux with neutron en- 
ergies up to a few GeV expected at 3200 mwc deep 
due to cosmic-ray /z interactions in the surround- 
ing rock and 30-cm lead shield and the resulting 
response of the CLOVER to the flux of 7 rays aris- 
ing from this flux. This simulation permits us to 
estimate rates in detectors situated in underground 
laboratories. 

The first two of these simulations are to verify the code's 
predictive power. The third is to aid in understanding the 
utility of depth to avoid neutron-induced backgrounds. 

The simulation package GEANT3-GCALOR [H, El 
is described in detail in Ref. [Tfj. In general, (n, n') 
reactions leave the target nucleus in a highly excited 
state which subsequently decays via a 7-ray cascade to 
the ground state. In the simulation, inelastic scattering 
cross sections for excitation to a given level depends on 
the properties of the ground and excited states. These 
cross sections were calculated using in-house-writtcn 
code based on Hauser-Feshbach [2(| theory modified by 
Moldauer [2l[. The validation of the Hauser-Feshbach 
theory has been the subject of several studies [22I I23I [2~i| . 
The simulated 7-ray flux arises from the relaxation of the 
initial excited-state distribution, which includes a large 
number of levels (60 states for 208 Pb(n, n'7) reactions, 
for example). The nuclear levels and their decay chan- 
nels were provided by the ENSDF[2f| database through 
the GEANT package. Note however, that the simula- 
tion did not predict every possible transition. In particu- 
lar the important 2041-keV and 3062-keV emissions from 
Pb were not part of this simulation. This situation arises 
because the simulation packages only have (n, n') cross 
sections for the lowest lying excited states for most nu- 
clei. It is set to zero for most other levels. The details of 
this simulation are described in detail in Ref. 10]. Here 
we study the effectiveness of the simulation to predict 
spectra resulting from (n,n'7). 

The simulation was done by generating neutrons with 
the appropriate energy spectrum outside the lead shield 
and propagating them through the shield including sec- 
ondary interactions that may add to the neutron flux and 



alter the energy spectrum. Fig. [3] shows a comparison 
between the data and the simulations for the CLOVER 
background run and AmBe run. Note only neutrons 
as primary particles were simulated for this comparison 
and the dominant difference between the two spectra is 
due to the room's natural radioactivity and non-neutron 
/i-induced processes. Here we excluded those processes 
from the simulation to emphasize the spectral shape, in- 
cluding lines, that are a direct result of neutron inter- 
actions. The similarity of the spectra in Fig. [3] indicates 
that the measured background spectrum is dominated by 
neutron-induced reactions. 

The uncertainty in the simulation is calculated by com- 
paring the well known peaks in Table |TT] which shows 
a comparison of the simulation to the line production 
for both background and AmBe runs. The measured 
neutron produced lines are within about 5% of the pre- 
dicted values from simulation, as is the continuum rate 
in the AmBe data. Therefore the (n,n'7) rates are well- 
simulated for nuclear states with well-defined cross sec- 
tions. (The continuum for the background data includes 
processes that were not simulated and hence is not a 
good measure of the uncertainty.) Because the neutron 
flux estimates come from these line strengths (See Sec- 
tion HIT]), the uncertainty in the flux cancels in these es- 
timates. The uncertainty in the measured neutron flux 
and spectrum underground(«35%) constrains the preci- 
sion to which such simulations can be verified and is well 
described in Ref. [l(|- This 35% uncertainty due to the 
flux is much larger than the uncertainty for the 7-ray 
line production described above. Therefore, a total un- 
certainty of 35% is used for all predictions of line rates 
underground throughout this paper. 



III. THE NEUTRON FLUX 

In this section, we use the data to determine the neu- 
tron fluxes we observed during our various experimental 
configurations. We then compare our measured cosmic- 
ray induced flux with that predicted from past measure- 
ments and our simulation. 



A. Ge(n, n'7) Analysis 

Spectral lines that indicate neutron interactions in nat- 
ural Ge detectors have been studied previously. See Ref- 
erences [1^, [27], HH [29||, for example. In particular, the 
sawtooth-shaded peaks due to 72 ' 74 Ge(n,n') at 693 keV 
and 596 keV respectively are clear indications of neu- 
trons and have been used to deduce neutron fluxes [3(ij |. 
Operating Gc detectors in a low-background configura- 
tion, these lines can be used to help interpret the back- 
ground components. Recent double-beta decay experi- 
ments 0, || have constructed their detectors from Ge en- 
riched in isotope 76. Although an appreciable amount of 
74 Ge remained (14%), 70 < 72 > 73 G e are depleted. For such 
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FIG. 3: Comparison of the measured and simulated AmBe 
spectra for the CLOVER detector surrounded by 10 cm lead 
and 15 cm of moderator. The upper plot shows the energy 
range between 10 - 3100 keV. The lower plot shows the range 
470 - 830 KeV where the most significant (n, n' 7 ) lines can 
be seen. The simulated AmBe neutron spectrum was nor- 
malized to the AmBe source strength for 6.13-h live-time. 
The measured total neutron flux in the background spectrum 
(see Section fill C[l was used to normalize the simulated back- 
ground spectrum. Note, only neutrons as primary particles 
were simulated for this comparison and the difference between 
the spectra is due to the room's natural radioactivity and non- 
neutron, //-induced processes. 



detectors, only lines originating in isotopes 74 and 76 
are useful for neutron interaction analysis. As these ex- 
periments reach for lower background, neutron-induced 
backgrounds become a greater concern and the diagnos- 
tic tools more important. 

Neutrons from (a,n) and fission reactions have an en- 
ergy spectrum with an average energy similar to the 
AmBe spectrum used in this study. Furthermore, the 
average energy of the AmBe neutrons is similar to that 
of the neutrons within the hadronic cosmic-ray flux im- 
pinging on our surface laboratory although the latter ex- 
tend to much higher energies. Therefore the Ge-detector 
signatures indicating the presence of neutrons described 
above will be similar to those arising from neutrons origi- 
nating from the rock walls of an underground laboratory. 
However, low-background experiments that use Ge de- 
tectors are typically deep underground and are shielded 



TABLE II: A comparison of the simulated to measured rates 
(Background: per hour and AmBe: Hz) for several lines pro- 
duced by neutron interactions. The 2041-keV and 3062-keV 
lines are not included in the simulation. 



Process 


7 -ray 


Backj 


jround-CLOVER 




Energy 


Simulation 


Measurement 


Y4 Ge(n,n' 7 ) 


596 keV 


56.21 




59.90(30) 


74 Ge(n,n' 7 ) 


254 keV 


2.63 




2.76(7) 


76 Ge(n,n' 7 ) 


2023 keV 


3.2x10" 


-7 


below sensitivity 


206 Pb(n,n' 7 ) 


537 keV 


4.82 




5.12(9) 


207 Pb(n,n' 7 ) 


898 keV 


6.21 




6.28(10) 


206 Pb(n,n' 7 ) 


1706 keV 


0.69 




0.72(3) 


206 Pb(n,n' 7 ) 


2041 keV 


none 




not seen 


Continuum region 


2000-2100 keV 


110.2 




187.35(19) 


207 Pb(n,n' 7 ) 


3062 keV 


none 




not seen 


Process 


7-ray 


AmBe-CLOVER 




Energy 


Simulation 


Measurement 


Y4 Ge(n,n' 7 ) 


596 keV 


1.8 




1.87 


74 Ge(n,n' 7 ) 


254 keV 


0.36 




0.41 


76 Ge(n,n' 7 ) 


2023 keV 


8.5x10" 


-4 


below sensitivity 


206 Pb(n,n' 7 ) 


537 keV 


0.15 




0.16 


207 Pb(n,n' 7 ) 


898 keV 


0.14 




0.20 


206 Pb(n,n' 7 ) 


1706 keV 


0.04 




0.04 


206 Pb(n,n' 7 ) 


2041 keV 


none 




not seen 


Continuum region 


2000-2100 keV 


7.01 




7.33 


207 Pb(n,n' 7 ) 


3062 keV 


none 




0.01 



from environmental radioactivity by a thick shield. This 
shield, typically made of Pb, is then usually surrounded 
by a neutron moderator. This configuration is effective at 
greatly reducing the neutron flux originating from (a,n) 
and fission reactions in the cavity walls of the under- 
ground laboratory. In contrast, although neutrons origi- 
nating from /i interactions underground are much rarer, 
they have much higher energy. Therefore these //-induced 
neutrons can penetrate the shield more readily and be- 
come a major fraction of the neutrons impinging on the 
detector. 



B. The AmBe Neutron Flux 

The estimate of the flux of neutrons with energies 
greater than 692 keV is given by [3(| HH, [HJ 



= k— , 

n y, 



(1) 



where I is the counts s" 1 under the asymmetric 692-keV 
peak, V is the volume of the detector in cm 3 (566 cm 3 ) 
and k is a parameter found by Ref. [3(| to be 900 ±150 
cm. For the 15-cm moderator data, this formula predicts 
a neutron flux of 2.3/(cm 2 s) whereas our simulation, 
using the known flux of the source, predicts 1.8/ (cm 2 
s). This difference (20-30%) is somewhat greater than 
the 17% uncertainty claimed by Ref. [3(J. The geometry 
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for our measurement was complicated and perhaps this 
added complexity of neutron propagation contributes to 
the difference. For the uncertainty associated with the 
flux of neutrons produced from cosmic ray [i, we use the 
35% value as it is much larger than the value associated 
with Eqn. [TJ 

For the Am-Be neutron source, the rate in the 692-keV 
peak is 2.406 ± 0.008 Hz. This results in $™ ihe = 3.8 
± 1.1 /(cm 2 s). This rate is an average over the two 
moderator configurations. The neutron flux during the 
10-cm moderator run is estimated to be about a factor 
2.3 larger than for the 15-cm moderator run. For the 
PopTop-AmBc run on Pb for the raw data (in coinci- 
dence with the Nal detector), the effective flux was 8.6 
± 2.6/(cm 2 s) (0.26 ± 0.08 /(cm 2 s)). 



C. Cosmic-ray Induced Neutron Fluxes 

In the background spectrum the rate in the 692-keV 
peak is 87.7 ± 0.4/hr. Using Eqn. ((TJ) with / = (2.44 ± 
0.06) x 10~ 2 Hz for the background spectrum, one obtains 
a fast neutron flux of § b * ck = (3.9 ± 1.2) xlO" 2 /(cm 2 
s) at the detector in our surface laboratory. 

Ref. [30j | provides a similar formula to estimate the 

thermal neutron flux, which is accurate to approximately 

30%. Using the intensity of the 139.68-keV 7-ray line of 
75m Ge . 



980/ 



139.68 



( e 139.68 



1.6)V" 



with 



fc 139.68 



1 - 



1 - e 



-V 1 ' 3 



(2) 



(3) 



where / = 47.2 ± 0.3 /h = 0.013 Hz is the event rate in 
the peak of 139.68-keV line and V is the volume of the 
detector in cm 3 . Using V = 566 cm 3 we obtain $^ c ' c = 
(9.1 ± 2.7) xl0~ 3 /(cm 2 s). We also measure the thermal 
neutron flux for the Am-Be neutron source, ^j™ 66 = 1.6 
± 0.5 /(cm 2 s). 

Thus the total neutron flux incident on the Ge detector 
measured for the background run is approximately 



back 



ifyback 



(4.8 ± 0.7) x 10" 2 /(cm 2 s). 



D. Neutron Flux as a Function of Depth 

In our basement laboratory, there are 3 primary 
sources of environmental neutrons. The largest contri- 
bution comes from the hadronic cosmic ray flux. The 
next largest arises from /i interactions in the 77 g/cm 2 
thick overhead concrete layer in the building. Finally 
there is the negligible contribution from (a,n) and fission 
neutrons from natural radioactivity in the room. The at- 
mospheric depth at the altitude of our laboratory is 792 
g/cm 2 . Including the concrete, the depth is 869 g/cm 2 . 



Using the analysis of Ziegler [33j, l34[ , the flux at our lab 
due to the hadronic flux can be estimated to be 3.0 times 
larger than that at sea level. The flux at sea level has 
been measured to be 1.22 x 10~ 2 /(cm 2 s) [HI result- 
ing in a flux in our laboratory of 3.7 x 10 -2 /(cm 2 s). 
To estimate the additional neutron flux originating from 
fi interactions in the concrete above our laboratory, we 
rely on our simulations of neutron generation and prop- 
agation. The simulation predicts 1.4 x 10" 2 /(cm 2 s) 
(3.3 x 10~ 2 /(cm 2 s)) for the muon-induced (hadronic) 
neutron flux inside the lead shield for a total simulated 
neutron flux of 4.7 x 10~ 2 /(cm 2 s) in acceptable agree- 
ment with our measurement of (4.8 ± 2.2) x 10~ 2 /(cm 2 
s) = (1.51 ± 0.69) x 10 6 ± /(cm 2 y). The success of this 
simulation lends credence to the neutron flux estimate in 
the following sections. 

The neutron flux onto the detector will be increased 
due to the neutron interactions with shield materials and 
neutron back-scattering from the cavity walls. For exam- 
ple, our simulations show that the fast neutron flux will 
increase by a factor of f=a 10 by traversing a 30-cm lead 
layer. Also, neutrons will backscatter from the cavity 
walls and reflect back toward the experimental appara- 
tus, effectively increasing the impinging neutron flux by 
a factor of 2-3 depending on the specific geometries of 
the detector and experimental hall. Therefore, it is im- 
portant to account for these effects when estimating the 
neutron flux at the detector. 

Muon-induced neutron production in different shield- 
ing materials and in the detector itself was also stud- 
ied in Ref. [Tcj. For example, with 30 cm of lead sur- 
rounding a CLOVER-style detector at 3200 mwe, the to- 
tal muon-induced neutron flux impinging on the detector 
was calculated to be (8.6 ± 4.0) x 10~ 8 /(cm 2 s) = 2.7 
± 1.2/ (cm 2 y). Some of the interactions resulting from 
these neutrons would be eliminated by a (i veto. Assum- 
ing a veto efficiency of 90% for muons traversing this lead 
shield, the effective neutron flux is estimated to be (2.0 
± 0.9) x 10~ 8 /(cm 2 s) = 0.63 ± 0.29 /(cm 2 y). The 
energy spectrum of neutrons at the lead/detector bound- 
ary at 3200 mwe is shown in Fig. [4] and has an average 
value of 45 McV. 

The average energy of the /x-induced neutrons is 100- 
200 McV and much higher than that of (a,n) neutrons 
(« 5 MeV). The simulated flux of the /i-induced neutrons 
((2 ± 0.9) x 10~ 8 /cm 2 s = 0.63 ± 0.29/cm 2 y) inside 
the detector shield at a depth of 3200 mwe is a factor 2.4 
greater than the simulated (a,n) flux surviving the shield 
((0.85 ± 0.39) x 10~ 8 /cm 2 s = 0.27 ± 0.12/cm 2 y). The 
average energy of these (a,n) originating neutrons is 3-5 
MeV at the detector surface. 

With this estimate of the neutron flux at depth and 
with our measurements of the Pb and Ge neutron- 
induced detector response, we can proceed to estimate 
these processes in underground Ge-detector experiments. 
There arc effects in addition to the incident flux, however, 
that must be taken into account when extrapolating our 
surface laboratory results to different geometries and lo- 
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FIG. 4: The effective neutron flux onto the simulated detector 
described in the text at a depth of 3200 mwe. Shown are the 
neutron flux from two sources: (1) the effective neutron flux 
induced by muons that transverse the surrounding rock and 
shielding materials assuming a 90% muon-veto efficiency and 
(2) the neutron flux from (a,n) reactions in the rock. 



cations. 

1. As the thickness of the Pb shield increases, addi- 
tional secondary neutrons will be generated. Our 
simulation predicts that a factor k s hi e id = 2.16 
more neutrons will be produced by a 30-cm thick 
shield as compared to a 10-cm thick shield. 

2. As the energy of the neutrons increases, the number 
of multiply scattered neutrons increases and there- 
fore the number of interactions that might produce 
a 7 ray increases. For the average energy of neu- 
trons at our surface laboratory (at 3200 mwe), the 
average scattering length is \l~7.1 cm (\ug = 12-5 
cm). 

3. Also as the energy increases, the number of states 
that can be excited in the target nucleus increases. 
In the shield at our surface laboratory (at 3200 
m.w.e), the average neutron energy is 6.5 MeV (45 
MeV). 

All of these factors can be incorporated into a scaling 
formula derived from our simulation. The rate (R^qi) 
of background near the region of interest (ROI) in an 
underground laboratory can be related to that measured 
in our surface laboratory (Rroi) as 
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~~ : ) K shield{{ 



E 



UG 



$1 



-R 
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(4) 

where $^ (^n°) ^ s the neutron flux in our surface labo- 
ratory (at 3200 mwe), E n is the neutron energy and E x 
is the excitation energy for a typical level. This formula 
reproduces our simulated results well and the uncertainty 
of its use is dominated by the precision of simulation. Us- 
ing the 2.6-MeV level in 208 Pb as an example, R%^>j ~ 
1.7 x 10- 5 R^ OI . Fig.EJshows a comparison between the 



Monte Carlo simulation and the scaling formula Eqn. [4] 
for several lines. 
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FIG. 5: The comparison between the Monte Carlo simulation 
of a detector as described in the text and the scaling formula 
for several excitation lines. The 35% uncertainties shown in 
this figure arise from the cross section uncertainty and the 
statistical uncertainty of determining the peak counts in the 
simulated spectra. The latter dominates. 



IV. ANALYSIS 

A. Ph(n, n'7) analysis 

If fast neutrons are present, then one will also see 7-ray 
lines from Pb(n, n') interactions. In very low background 
configurations, 7 rays from neutron-induced excitations 
in 208 Pb and 207 Pb can be masked by or confused for de- 
cays of 208 T1 and 207 Bi respectively. Therefore it is the 
stronger transitions in 206 Pb (537.5, 1704.5 keV) that are 
most useful for determining if these processes are taking 
place. In 207 Pb the relative strength of the 898-keV tran- 
sition, compared to the 570- and 1064-keV transitions, is 
much stronger when it originates from 207 Pb(n,n'7) as 
opposed to 207 Bi (3 decay to 207 Pb. Therefore this line 
can also be used as a tell-tale signature of neutron inter- 
actions. 

Our data show indications of 206,207 ' 208 Pb(n, n'7). As 
noted earlier, the 2614-keV 7 ray from 208 Pb can orig- 
inate from 208 T1 decay or from 208 Pb(n,n'7). The 692- 
keV peak arises only from neutron interactions on 72 Ge. 
Since the Pb shielding was similar in both the background 
and AmBc runs, we can compare the ratio of the rate in 
the 2614-keV peak to that in the 72 Ge(n,n'7) 692-keV 
peak in the two data sets to deduce the fraction of the 
2614-keV in the background run that can be attributed 
to neutron interactions. This ratio in the Am-Be spec- 
trum is 0.30 and that in the background spectrum is 0.45, 
and therefore, we conclude that f=a 67% of the strength 
in the background run is due to neutron reactions and 
the remainder is due to 208 Tl decay. Clearly, in our sur- 
face laboratory, environmental neutrons are a significant 
contributor to the observed signal. 



11 



B. The Special Cases of the 2023-keV, 2041-keV 
and 3062-keV 7 rays 

The 2023-keV level in 76 Ge can be excited by neutrons. 
The simulation predicts that this line is too weak to be 
seen in the CLOVER AmBe data, but the CLOVER is 
built of natural Ge. In the enriched detectors planned 
by future double-beta decay experiments, the fraction of 
isotope 76 is much larger and this line would be enhanced. 
Still our simulation (Tabic [Vj predicts it would be a very 
small peak. 

The 3714-kcV level in 206 Pb can emit a 2041-keV 7 
ray. We only observed a candidate 7-ray peak in the 
coincidence data (Ge detector event in coincidence with 
a 4.4-McV 7 ray in the Nal detector) with the AmBe 
source radiating the Pb shield around the PopTop detec- 
tor. The magnitude of this peak if it exists is small and 
not convincing. We use this data set to place a limit on 
the production rate of this line as it results in the most 
conservative limit. 

In the AmBc-irradiated CLOVER and the non- 
coincidence PopTop spectra, we observed a 3062-keV 7 
ray that we assign to a transition from the 3633-keV level 
in 207 Pb. (See Fig. [6]) This line is only present when 
Pb surrounds the detector: It is absent when Cu forms 
the shield. The statistical sensitivity was too weak in 
the PopTop coincidence spectrum to observe this weak 
line. In the CLOVER AmBe spectrum, the rate of this 
line is 5.3 x 10~ 3 that of 596-keV 74 Ge peak rate and 
in the raw AmBe PopTop with Pb spectrum the ratio is 
4.5 x 10~ 3 . From these data we can estimate an approx- 
imate rate that these dangerous backgrounds would be 
produced for a given neutron flux. In our surface labora- 
tory, the CLOVER background rate for the 74 Ge 596-keV 
peak was 59.9 events/hr. This leads to a predicted rate of 
0.3 events/hr in the 3062-keV peak. Note that our data 
indicate that any peak at 3062 keV is statistically weak 
(< 0.2 events/hr) but reasonably consistent with this pre- 
diction. Note, in the AmBc-CLOVER runs, polyethylene 
blocks were used to increase the flux of thermal neutrons. 
It appears that these blocks contain some CI and there- 
fore we see indications of Cl(n,7) lines. Even though 35 CI 
has a neutron capture line at 3062-keV, we do not assign 
the observed line in our data to that process. Because the 
35 Cl(n,7) line at 2863 keV is not observed and because 
the line at 1959 keV is weak, we conclude that assigning 
this line to CI would be inconsistent with the predicted 
line ratios for neutron capture. However, a concern re- 
garding our assignment of the 3062-keV line to 207 Pb is 
that one also expects a 2737-keV emission from the same 
3633-keV level. This companion 7 ray is not observed 
in our data and we plan future measurements dedicated 
to measuring the neutron-induced relative intensities of 
these two lines. If we assume that the entire rate (0.01 
Hz) of the 3062-keV line is due to (n, n'7), we can make a 
crude estimate of the cross section by scaling to the rate 
in the 2614-keV line. The cross section for the 2614-keV 
(n.n'7) is 2.1 b ± 10% [23|. Using this cross section, the 



relative rates in the two peaks, the different isotopic ra- 
tios of 208 Pb and 207 Pb, and the different 7-ray detection 
efficiencies, the average cross section for 207 Pb(4.5-MeV 
n, n' 3062-keV 7 ray) is estimated to be 75 mb. The 
uncertainty is estimated to be about 20% or ~ 15 mb. 




3080 3100 
Energy (keV) 



FIG. 6: The energy spectrum near 3000 keV showing the 
207 Pb(n, n'7) 3062-keV 7-ray line in the AmBe spectrum with 
the CLOVER. 



From measurements with the CLOVER and a 56 Co 
source, which has 7-ray energies near 3100 keV, we expect 
0.13 DEP events per full-energy 7-ray event. Therefore, 
in our surface lab, we expect 0.03 events/hr in the DEP 
at 2039 keV due to 207 Pb(n, n'7). This is well below 
our continuum rate of 2.5 events/keV-hr or 10/hr in an 
energy window corresponding to a 4-keV wide peak. 

Since our simulation does not predict all these lines, we 
summarize the measured rates normalized to the neutron 
flux in Table HLT1 to provide simple scaling to different ex- 
perimental configurations. The uncertainties in Table ITTT1 
are estimates based on a minor contribution of the sta- 
tistical uncertainty in the peak strengths and a major 
contribution resulting from the 35% uncertainty in the 
neutron flux determination as described in Section IIII Bl 
Because the uncertainty is mostly systematic, there is a 
good possibility that the total uncertainties for each in- 
dividual measurement are correlated. Therefore, to esti- 
mate the average values in this Table, we took a straight 
average of the individual values and then assigned an un- 
certainty equal to the largest fractional value. This pro- 
cedure, although not rigorous, is more conservative than 
a weighted average. In addition, some peaks were not 
observed in all spectra. The upper limits on the strength 
of these peaks were estimated from the rates of weak- 
est peaks observed near the associated energy region in 
the spectrum. Such peaks are considered to represent 
the level of sensitivity of our peak detection procedures. 
The 2041-keV line is a special case. We quote an up- 
per limit based on the only spectrum that indicated a 
possible peak. 

These measurements were done for a CLOVER de- 
tector inside a 10-cm Pb shield. The relative energy- 
dependent efficiency (e r ei) for a full-energy peak in the 
CLOVER can be approximated by, 
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TABLE III: The raw count rates for select processes normalized to neutron flux of l/(cm 2 y). When extrapolating to neutron 
energies distant from that near the measurements, the uncertainty (35%) associated with the extrapolation must be included. 
See text for a discussion of the uncertainty estimates in this table especially with respect to the average. 



Process Rate Rate Rate Rate 

AmBe-CLOVER AmBe-PopTop Background Average 

/ c vents \ 



2Ub Pb(n,n'537- keVf) 


13.9±4.9 


20.5±7.2 


12.2±4.3 


15.5±5.4 


74 Ge(n, n'596 - keV-f) 


164±57 


unresolved" 


142±50 


153±54 


207 Pb(n,n'898 - keVj) 


17.5±6.1 


21.2±7.4 


14.9±5.2 


17.9±6.3 


206 Pb(r 1 ,n'1705 - keVj) 


3.6±1.3 


3.2±1.1 


1.7±0.6 


2.8±1.0 


206 Pb(n,n'2041 - keV-y) 


<6.3 


< 0.5 


<1.7 


< 5.0 6 


207 Pb(ri,n' 3062 - keVy) 


0.88±0.31 


0.6±0.21 


<1.0 


0.7±0.2 


Continuum Rate from Pb, Ge(n, n'7) 


2.6±0.9 


( events \ 
V (keVty) > 

2.0±0.7 


2.5±0.9 


2.4±0.8 



"In the PopTop data, the 596-keV line was not resolved from 
nearby lines. 

'The 2041 line was not observed in any of our spectra, however, 
a weak peak-like feature was present in the AmBe-PopTop coinci- 
dence data. We used the upper limit for the rate in that peak as 
the "average" as we considered this to be most conservative. 



-(E„-148) 

e re i = 0.15 + 0.93e rm — , (5) 

where E~ is the 7-ray energy in keV. This expression is 
normalized to 1.0 at 209 keV and is estimated to have an 
accuracy of about 20% near 200 keV improving to about 
10% at 2600 keV. The quoted relative efficiency for each 
of the 4 individual CLOVER detectors is 26% at 1.33 
MeV as quoted by the manufacturer. Table [LTTl does not 
incorporate this efficiency correction, therefore the table 
presents the measured count rates with a minimum of 
assumptions. The thickness of Pb is large compared to 
the mean free path of the 7 rays of interest, therefore, the 
scaling should hold for other thick-shield configurations. 
Even so, the rates will be geometry-dependent so these 
results can only be considered guides when applied to 
other experimental designs. The rate of these excitations 
also depends on neutron energy. For the background run 
(AmBe run) the average neutron energy is « 6.5 MeV (« 
5.5 MeV). Our simulations predict that the rate of these 
excitations scales as energy to the 0.81 power. 

V. DISCUSSION 

A. A Model of the CLOVER Background 

We can use these experimental results to create a back- 
ground model for our surface lab and deduce the contri- 
bution to the continuum near 2039 keV due to (n,n'7) re- 
actions. We then use simulation of high-energy neutron 
production and propagation to extrapolate this model 
to better understand experiments done at depth. The 
measured rate for the continuum near 2039 keV was 14.8 



events/ (keV kg d). For the Th-wire data, this continuum 
rate was 0.10 events/ (keV s) (2900 events/ (keV kg d)) 
and for the AmBe data it was 0.09 events/(keV s) (2600 
events/ (keV kg d)). To determine the neutron-induced 
continuum rates in the AmBe data, however, we have 
to correct for the contribution from the tail of two high- 
energy 7 rays that are not part of the neutron-induced 
spectrum in the background. These are the 7 rays from 
the 2223-keV p(n, 7)d and the 4.4-MeV 7 rays originat- 
ing from the (a,n) reaction of the AmBe source itself. 
Although only w 10% of these AmBe 7 rays penetrate 
the 5-cm Pb shield, there is still a significant flux. 

A simple simulation of the detector response to 2223- 
keV 7 rays can easily determine ratio of the rate in the 
2039-keV region to that in the full-energy peak. Simula- 
tion indicates that this ratio is 5.2 x 10~ 2 /keV. Since 
the full-energy peak count rate is 5.813 Hz, we find this 
contribution to the continuum is 0.03 events/(keV s). For 
the 4.4-McV 7 ray from the source itself, simulation must 
determine an absolute rate in the continuum because the 
high-energy threshold prevented the observation of the 
full-energy peak or its escape peaks. The simulation pre- 
dicts 0.03 events/(keV s). Subtracting these two contri- 
butions from the continuum rate for the AmBe source 
near 2039 keV results in a final value of 0.03 events/ (keV 
s) or 860 events / (keV kg d) . 

Our background measurements were done without a 
cosmic ray anti-coincidence system. From auxiliary mea- 
surements with a scintillator in coincidence with the 
CLOVER and a similar shielding geometry, we measured 
the rate of \x passing through the detector. In the contin- 
uum near the 2039-keV region, the rate is 5.4 events/ (keV 
kg d). 

From the Th-wire source data, we measure the ratio 
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TABLE IV: A summary of the count rate in the CLOVER 
background data in the energy region near 2039 keV based 
on the model deduced for the surface lab described in the 
text. The precision of the neutron-induced and muon-induced 
spectra simulations fSection lll Cl and Ref [l(|) is estimated to 
be about 35%. We take this to be a conservative estimate for 
the uncertainties associated with this Table. 



Process 


CLOVER Event Rate 




Surface Lab 




events/(keV kg d) 


neutron-induced 


8.3±2.9 


208 Tl Compton scattering 


0.7±0.3 


high energy /j, continuum 


5.4±1.9 


Total from model 


14.4±5.0 


Measured Rate 


14.8±0.2 



of the rate in the continuum near 2039 keV to that in 
the 2614 keV (16.3 Hz) peak to be 6 x lCT 3 /keV. Of 
the 2614-keV peak rate in the background data, pa 33% 
is due to 208 T1 decay. Scaling from the 2614-keV peak in 
the background data, the count rate near the 2039-kcV 
region due to the Compton tail of the 208 T1 2614-keV 
peak is 0.7 events/(keV kg d). 

The remainder of the 2614-keV peak is due to neutron- 
induced processes. The contribution due to neutrons can 
be estimated from the AmBe data. For the AmBe data, 
the ratio of the rate in the continuum near the 2039- 
keV region (0.03 events(keV s)) to that for the 596-keV 
74 Gc(n,n'7) peak (1.87 Hz) is 1.6 x 10~ 2 /kcV. Scaling 
from the 74 Ge peak rate in the background data (59.9/h) 
indicates a rate of 7.8 events/ (keV kg d) in the contin- 
uum near 2039 keV. That is, 53% of the events in that 
region are due to neutrons. One can do a similar scal- 
ing from the 692-keV 72 Ge rates. Here the ratio is 1.3 x 
10~ 2 /keV and the continuum rate is 8.8 events/(keV kg 
d). We use the average of the Ge values as our estimate 
(8.3 cvcnts/(kcV kg d) = 3030 events/(kcV kg y)) for 
the neutron induced contribution to the continuum rate. 
Table IIVI summarizes the deduced contributions to the 
spectrum in the 2039-keV region in the CLOVER back- 
ground spectrum and the following section discusses how 
these data are used along with simulation to estimate 
rates in experimental apparatus underground. 



B. Solving the Problem with Overburden 

The primary purpose of this study is to better under- 
stand the impact of neutrons on the background for fu- 
ture double-beta decay experiments. In this subsection, 
we use neutron fluxes from our simulations of the surface 
laboratory, measurements with the AmBe source, and 
simulations of the neutron flux in an underground lab- 
oratory to estimate the contribution of neutron-induced 
backgrounds underground. In the following subsection, 



we examine data from previous underground experi- 
ments. 

The simulation of neutron processes in the 10-cm 
Pb shield and Ge comprising the CLOVER detector at 
the altitude of our laboratory predicts about 1594±558 
events/ (keV kg y) between 2000 and 2100 keV due to lead 
excitation and about 1337±468 events/ (keV kg y) in this 
energy region due to germanium excitation. Our mea- 
sured value for the neutron-induced events is 3030±1061 
events/ (keV kg y) to be compared with this predicted 
value of 2931±1026 events/(keV kg y). 

The simulation of the CLOVER within a 30-cm 
lead shield at a depth of 3200 mwe, predicts about 
0.019±0.007 events/(keV kg y) contributed from lead 
excitation and about 0.016±0.006 events/(keV kg y) 
contributed from germanium excitation for a total of 
0.035±0.012 cvents/(kcV kg y). One can also just 
scale our surface-laboratory measurement of the neutron- 
induced rate near 2039 keV by the factor derived from 
Eqn. H] above. This results in 0.05±0.02 events/(keV 
kg y). For a detector like the CLOVER, analysis based 
on pulse shape discrimination (PSD), and the response 
of individual segments or crystals can help reduce back- 
ground based on its multiple-site energy deposit nature. 
These backgrounds can then be distinguished from the 
single-site energy deposit character of double-beta decay. 
We have measured the background reduction factor via 
these techniques to be 5.9 for the CLOVER [ljj. 

Reference fioj provides a quick reference formula to 
estimate the neutron flux as a function of depth. The 
fi flux and its associated activity is reduced by ps 10 for 
each 1500 m.w.e of added depth. Future double-beta 
decay experiments hope to reach backgrounds near 0.25 
events/ (keV t y). Our estimate of the rate at 3200 mwe 
is 35-50 events/ (keV t y), which is a factor of 150 above 
the goal. Hence, greater depths would be desirable. 



C. Discussion of Previous Underground 
Experiments 

Previous Gc-based double-beta decay experiments con- 
ducted deep underground [f| @ set the standard for low 
levels of background. The future proposals however [1, Q 
hope to build experiments with much lower backgrounds. 
In this subsection, we estimate how large the neutron 
contribution was to the previous efforts and future de- 
signs. Using the scaling summarized in Eqn. 01 we can 
compare the expectations of our simulated underground 
apparatus with previously published results. Table IVl 
shows this comparison. This table also presents a sum- 
mary of how the rates would be affected by a change in 
depth only. The IGEX collaboration Q has not pub- 
lished its data in sufficient detail to do a similar com- 
parison. Other underground Ge detector experiments do 
not have the required sensitivity. 

The Heidelberg-Moscow experiment @ is a critical 
case study for such backgrounds and it was operated at 
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TABLE V: A summary of the key count rates arising from neutron interactions in the CLOVER background data in the energy 
region near 2039 keV as predicted by our analysis for three representative depths. The shield thickness is taken to be 30 cm 
and a veto system with an assumed efficiency of 90% is included. Except for the 2023-keV line, we used the scaling of Eqn.|4]to 
scale our CLOVER background measurements to the 3200 mwe depth and then used the muon fluxes at WIPP, Gran Sasso, and 
SNOLAB[i(| to scale to the other depths. The Ge rates are also scaled for an enriched detector (86% isotope 76, 14% isotope 
74). The scalings require the results from the simulations. The uncertainty is dominated by the simulated flux uncertainty 
and is estimated to be 35%. Since we did not observe the 2023-keV line, we used simulation to predict the rate. We used the 
measured upper limit for the 2041-keV line. For comparison, the results of Ref. [13] is shown. Reference [13] claims a result for 
zero-neutrino double-beta decay in an experiment performed at 3200 mwe. We entered the claimed event rate for that process 
in the same row as the 2041-keV line for comparison. The rate limits for the other lines assigned to Ref. [37]] result from our 
estimates based on the figures in their papers and does not come directly from their papers. 



Process 



1600 mwe 



3200 mwe 



6000 mwe 



74 Ge 596 keV 
76 Ge 2023 keV 
206 Pb 537 keV 

207 



Ref. [37] 



Pb 898 keV 
3 Pb 1705 keV 
Pb 2041 keV 
207 Pb 3062 keV 
continuum 



2()l> 



19400±6800/(t y) 

30±10/(t y) 
4400±1500/(t y) 
5300±1900/(t y) 

610±210/(t y) 
<1300±450/(t y) 

145±51/(t y) 
880±310/(keV t y) 



1130±400/(t y) 
2±l/(t y) 
250±88/(t y) 
310±110/(t y) 
36±13/(t y) 
<74±26/(t y) 
8.4±2.9/(t y) 
50±17/(keV t y) 



15±5/(t y) 
0.02±0.01/(t y) 
3.4±1.2/(t y) 
4.2±1.5/(t y) 
0.5±0.2/(t y) 
<1.0±0.3/(t y) 
0.1±0.03/(t y) 



<800/(t y) 
300/(t y) 



400/(t y) 
<71/(t y) 



0.7±0.24/(keV t y) 110/(keV t y) 



a depth of 3200 mwe. One is clearly led to consider if 
the 3062-keV 7 ray can explain the signal reported in 
Ref. [H, HtJ ■ Figure 36 in Ref. [37j , shows that no more 
than a few counts can be assigned to a 3062-keV 7 ray. If 
the 23 counts assigned to double-beta decay were actually 
a DEP from this 7 ray, the one would expect 175 counts 
or so in the 3062-keV peak. Therefore, it is difficult to 
explain the claimed peak by this mechanism. It is also 
clear that the predicted rate of the 2041-keV 7 ray is too 
low to explain their data. The data from Ref. [37| show 
lines at 570 and 1064 keV and the authors assigned these 
lines to 207 Bi present in the Cu. However, the spectra dis- 
played in Fig. 13 of that paper shows that only detectors 
surrounded by Pb indicate the 570-keV line. Since there 
is no evidence for the 898-keV line in the data, we agree 
with the 207 Bi assignment, however, we hypothesize that 
it must reside in the Pb and not the Cu. Perhaps this 
contamination is cosmogenically produced in Pb when it 
resides on the surface and not as a result of bomb testing 
as hypothesized by the authors [37 1. 

Reference [13 also observed lines at 2011, 2017, 2022, 
and 2053 keV. These lines had rates of approximately 
500/ (t y), 500/ (t y), 300/ (t y) and 380/ (t y) respectively. 
The line at 2022 keV is near a line we predict at 2023 keV. 
Reference [3?| attributes these lines to weak transitions 
in 214 Bi. From our analysis it is indicated that a negligi- 
ble fraction of the peak at 2023 keV is neutron-induced. 
However, since the predicted strength of the tell-tale lines 
that would indicate a presence of neutron interactions is 
just below the sensitivity of that experiment, this conclu- 
sion is not without uncertainty. It has been pointed out 
that the strength in the 2022-keV line is too strong with 
respect to the 214 Bi branching ratios even when sum- 
ming uncertainties are taken into account [HI, [39| . The 
analysis in Ref. [38], however, was based on a incorrectly 



normalized Fig. 1 in Ref [6J. A recent analysis [39| tak- 
ing this into account still points to an inconsistency in 
the line strengths. This discrepancy could be resolved if 
one attributes a significant fraction of that peak to neu- 
tron interactions on 76 Ge. Such an attribution is not 
supported by our simulations. 

Reference I36| simulates the background in the 



Heidelberg-Moscow experiment resulting in a predicted 
signal of 646 ± 93 counts in the region between 2000 and 
2100 keV during an exposure of 49.6 kg-y. This is a count 
rate of 130/(kcV t y) to be compared with the quoted 
measured value for the data period simulated of 160/ (keV 
t y). Their estimate indicates that only 0.2/ (keV t y) are 
due to neutrons and they argue that /^-generated neu- 
trons are a negligible contribution. Our estimates indi- 
cate that neutrons are a more significant contribution and 
that the (i contribution is significant. We are aware of no 
direct neutron flux measurements for neutrons above 25 
McV. The flux of neutrons with energy greater than 25 
MeV is estimated in Ref. [H] to be 10~ n /(cm 2 s ) and 
they considered these neutrons to produce a negligible 
contribution to the background. In contrast, the simula- 
tion in Ref. gives 56 x 10- n /(cm 2 s) at 3200 m.w.e 
for the neutrons with energy greater than 25 MeV. We 
use this higher flux value and as a result, our estimate 
of the background rate near 2 McV of 50/(keV t y) is 
comparable to the excess (30/(keV t y)) of the measured 
rate in comparison to the simulated rate in Ref. [36| . 



D. Is Copper an Alternative to Lead? 

One has to consider the existence of a DEP line at the 
double-beta decay endpoint a serious design considera- 
tion for Ge-detcctor experiments. From the above anal- 
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ysis, the dangerous lines at 2041 and 3062 keV due to 
Pb(n, n'7) are not significant contributors to the spec- 
trum of Ref. [33] • However, as future efforts reduce 
the natural activity irradiating the detectors, these Pb- 
neutron interactions will become important. One solu- 
tion could be the use of Cu as a shield instead of Pb. 
Copper is rather expensive and building the entire shield 
of this material is probably not necessary. A thick inner 
liner of Cu might suffice, but if a peak is observed and 
Pb is present near the detector, arguments based on the 
spectrum near 3062 keV will be critical. 

Although the problematic lines we observed in the Pb 
data were absent in our Cu data, the shields were too 
dissimilar to make a quantitative comparison regarding 
the effectiveness of reducing the continuum background. 
Furthermore, our experience with the lead and the sim- 
ulation of (n, n'~f) spectra reduces confidence in the con- 
clusion regarding the Cu in the absence of such data. We 
are preparing better experimental studies to address this 
question. 

VI. CONCLUSION 

As double-beta decay experiments become more sen- 
sitive, the potential background must be constrained to 
ever-lower levels. Much progress has been made in reduc- 
ing naturally-occurring radioactive isotopes from mate- 
rials from which the detector is constructed. As these 
isotopes that have traditionally limited the experimen- 
tal sensitivity are eliminated, rarer processes will become 



the dominant contributors. Here we have considered 
neutron-induced processes and have quantified them. Re- 
actions involving neutrons can result in a wide variety of 
contributions to the background. That is, no single com- 
ponent is likely to dominate. Therefore, tell-tale signa- 
tures for neutrons are needed and were identified in this 
work. 

In addition to the general continuum background 
that neutrons might produce, two specific dangerous 
Pb(n, n'7) lines were identified. These two backgrounds 
can be significantly reduced using depth and/or an inner 
layer of Cu within the shield. In particular, the 3062- 
keV transition in 207 Pb has a double escape peak at the 
endpoint energy for double-beta decay in 76 Ge. A com- 
parison of past double-beta decay data indicates the rate 
of this transition is too small to explain a claim of double- 
beta decay. 
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